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ABSTRACT 

The  article  describes  the  modernization  of  Unmanned  Aerial  Vehicle  (UAV)  to  perform  obserrations  of  satellites.  The 
technical  requirements  applied  to  the  device  and  the  parameters  of  the  payload  and  its  parameters  are  verified.  Various 
models  used  for  aerodynamic  calculations  are  displayed  and  their  pros  and  cons  are  revealed.  Various  methods  of 
constructing  solid-state  models  are  also  displayed.  Several  different  shapes  of  the  dome  for  equipment  placement  are 
proposed.  The  optimal  shape  of  the  dome  was  chosen  based  on  the  results  of  calculations  on  the  pressure  distribution  on 
the  body.  The  angle  of  attack  was  chosen  based  on  the  calculation  of  the  resistance  and  lift  forces. 
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1.  INTRODU CTION 

According  to  the  military  doctrine  of  the  Republic  of  Kazakhstan,  modern  military  conAicts  have  been 
characterized  by  their  activity,  the  brevity  of  the  scale,  spheres  of  warfare  (space  and  information  space)  as  well  as 
high  voltage  strength  and  resources  of  the  state  in  the  armed  struggle  (Ministry  of  Defense  of  the  Republic  of 
Kazakhstan  official  resource).  In  this  regard,  it  is  of  great  importance  for  the  defense  of  the  country  to  get  the 
direction  associated  with  continuous  spacecraft  tracking  and  compiling  information  about  them.  Satellite  radar  is 
hampered  by  its  small  size  and  long  distances.  With  this,  we  put  forward  the  idea  of  monitoring  these  devices  using 
photomultipliers,  which  allow  us  to  study  even  faint  objects.  However,  observations  in  the  visible  range  from  the 
earth's  surface  are  often  hampered  by  dense  tropospheric  clouds.  Based  on  our  analysis  (Savinkin  et  al .,  2018),  it 
was  revealed  that  statistically  dense  clouds  are  formed  at  altitudes  not  exceeding  3.5  km. 

Such  heights  are  available  to  a  wide  range  of  aircraft  and,  in  particular,  unmanned  aerial  vehicles.  Based 
on  this,  we  see  a  solution  to  the  problem  of  improving  the  efficiency  of  monitoring  the  upper  sphere  in  the  design 
(and  possibly  adaptation)  of  an  unmanned  aerial  vehicle  capable  of  rising  to  heights  of  up  to  4  000  m  (4.0  km)  with 
the  possibility  of  installing  a  container  with  scientific  equipment.  This  container  can  be  a  dome  camera  with  a  gyro 
platform  for  panoramic  shooting  during  flight. 
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The  study  of  quite  a  large  number  of  different  aircrafts  led  to  the  choice  of  the  basic  model  of  the  UAV  GeoScan 
201,  the  characteristics  of  which  are  given  below  (Savinkin  et  al.,  2018): 

Max.  takeoff  weight:  8.5  kg 

Flight  speed:  64-130  km/h 

Flight  duration:  up  to  3  hours 

Max.  pay  load  weight:  1.5  kg 

Max.  flight  altitude:  4000  m  (4.0  km) 

Operating  temperature:  -  40  °C  to  +  40  °C 

As  can  be  seen  from  these  characteristics,  this  device  has  a  sufficient  practical  ceiling  to  overcome  the  dense 
tropospheric  clouds  and  is  also  able  to  take  quite  a  decent  mass  of  payload,  which  allows  us  to  modify  this  UAV  for  our 
tasks. 

As  a  payload,  we  selected  two  video  tracking  devices  -  one  camera  with  a  viewing  angle  of  1 80°  for  viewing  the 
front  hemisphere  of  the  aircraft,  installed  in  the  front  part  of  the  body,  and  the  second  aimed  at  the  zenith  and  equipped 
with  a  photomultiplier  for  tracking  satellites  installed  in  the  middle  part  of  the  fuselage. 

To  justify  the  choice  of  the  optimal  aerodynamic  shape  of  the  dome,  it  is  necessary  to  calculate  the  different  types 
of  dome  on  the  base  glider.  For  this  purpose,  various  mathematical  models  describing  the  flow  of  subsonic  air  masses  were 
chosen. 


The  Navier-Stokes  equations  are  the  basis  for  solving  aerodynamics  problems  (Belotserkovsky  et  al.,  2000; 
Belotserkovsky  et  al.,  2001;  Shevelev  and  Maksimov,  2000;  Shevelev  and  Maksimov,  2003a;  Shevelev  and  Maksimov, 
2003b).  However,  solving  the  complete  Navier-Stokes  equations  remains  a  very  difficult  task  even  for  geometrically 
simple  cases.  Solving  problems  of  external  aerodynamics  for  real  aircraft  can  be  done  only  with  the  use  of  specific 
simplifications.  Previous  researchers  -  Belotserkovsky  O.  M.,  Andrushchenko  V.  A.  and  Shevelev  Yu.  D.  used  three 
different  methods  of  solving  the  Navier-Stokes  equations  as  a  mathematical  model  for  the  study  of  viscous  incompressible 
gas  flow  near  bodies  of  complex  shapes:  (1)  the  method  of  "artificial"  compressibility,  (2)  the  method  of  splitting  by 
physical  variables,  generalized  to  the  case  of  arbitrary  curvilinear  coordinates  and  (3)  the  method  of  calculating  the 
Navier-Stokes  equations  in  the  “thin”  layer  approximation.  However,  in  the  proposed  mathematical  model, 
Belotserkovsky,  O.  M.,  Andrushchenko,  V.  A.  and  Shevelev,  Yu.  D.’s  emphasis  is  placed  on  the  principles  of  distribution 
of  velocity  Vx,  Vy,  Vz,  and  the  pressure  of  the  flow  medium  impact  px,  py,  pz  does  not  take  into  account  the  dynamic 
impact  of  body  parts  and  feature  of  the  layout  of  equipment  on  the  aircraft.  The  expansion  of  the  functionality  of  the 
developed  aircraft,  such  as  the  ability  to  observe  the  upper  hemisphere  of  the  outer  space  as  well  as  the  tracking  of  natural 
objects,  inevitably  entails  a  structural  change  that  affects  the  dynamic  loads,  namely  the  air  resistance  forces  Fr  acting  on 
the  apparatus.  Therefore,  there  is  a  need  to  improve  the  mathematical  model  of  the  study  of  the  distribution  of  forces  and 
pressure,  taking  into  account  the  dynamic  and  structural  changes  of  the  apparatus.  Also,  the  task  of  mathematical  modeling 
is  complicated  by  the  fact  that  a  significant  role  in  the  adequacy  of  the  model  is  played  not  only  by  established  processes, 
but  also  by  transients  in  the  boundary  layers  of  interaction  of  the  contact  medium  with  the  object  of  study,  for  example,  in 
the  approximation  of  a  thin  layer. 
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To  study  the  pitch  of  an  unmanned  aerial  vehicle,  it  is  important  to  describe  mathematically  the  process  of  flows 
of  yiscous  compressible  gas  near  bodies  of  non-regular  shape.  For  this  purpose,  various  methods  of  calculation  of 
stationary  and  non-stationary,  three-dimensional  Euler  and  Navier-Stokes  equations  written  in  conservative  form  are  used. 
To  solve  the  scientific  and  technical  problem,  several  solutions  are  proposed  based  on  finite  difference  and  finite  volume 
methods,  using  explicit  and  implicit  schemes  of  increased  order  of  accuracy  in  time  and  space,  such  as  TVD,  ENO,  etc. 

The  system  of  Navier-Stokes  equations  is  used  as  a  basis  for  improving  the  mathematical  model  of  the 
aerodynamic  system. 


2.  METHODS 

2.1  System  of  Nayier-Stokes  Equations 


Unsteady  Navier-Stokes  equations  for  three-dimensional  flow  of  compressible  gas  in  dimensionless  form  in  Cartesian 
coordinate  system  X  =  (x,  y,  z)  have  the  form 


^  +  £(f-^v)+^(F-Fv)+£(G-Gv)  =  0, 
where  functions 


p  PU 

pu  pu 2  +  p 

U  =  {pv},  E  =  {  puv  },  F 
pw  puw 

e  (e  +  p)u 


p\  pw 

puv  puw 

{pv2  +p},G  =  {  Pwv  }. 

pvw  pw2  +  p 

{e  +  p)v  ( e  +  p)w 


Here,  t  is  time,  [rho]  is  the  density,  (u,  v,  w)  -  components  of  the  velocity  vector  V,  respectively,  in  directions  (x, 
y,  z),  p  -  pressure,  e  -  total  energy  per  unit  volume  of  gas,  which  for  an  ideal  gas  can  be  represented  in  the  form 


(  2,2.  2\ 
U  +  V  +  w 
e  =  p  £  + - 


V 
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1  p 

— - -  the  internal  energy  of  a  gas,  [gamma]  -  the  ratio  of  heat  capacities 

r-ip 


at 


constant  pressure  and  volume. 
Further, 
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Dimensionless  variables  are  defined  through  dimensional  variables  that  are  marked  with  a  stroke: 

t=  ^x  =  x-,v  =  ^p  =  ^„P  =  Ar  =  -+  =  ^. 
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The  subscript  "zero"  means  the  value  of  the  parameter  in  the  unperturbed  stream.  Here,  L  is  the  characteristic  size, 
X  =  (x,  y,  z)  and  V  =  (u,  v,  w). 

It  is  assumed  that  the  Prandtl  number  Pr  =  is  constant,  here  cp  -  is  the  heat  capacity  coefficient,  [lambda]  -  is 
the  heat  conductivity  coefficient,  the  viscosity  coefficient  [mu]  can  be  calculated  by  the  formula  [i  =  Tk,  where  0.5  <  k  < 
1.5,  Re  =  yJp^oPo  Lf/iif0-  is  the  Reynolds  number. 

The  system  of  differential  equations  is  supplemented  by  the  equation  of  state  p  =  pRT,  where  T  -  is  the 
temperature,  R  -  is  universal  gas  constant.  In  dimensionless  form,  the  equation  of  state  has  the  form  p  =  pT.  Thus,  a 
complete  system  of  Navier-Stokes  equations  for  perfect  gas  flows  in  the  absence  of  external  mass  forces  is  given.  It  is 
assumed  that  the  heat  to  the  medium  can  only  come  as  a  result  of  thermal  conductivity. 

2.2  System  of  Navier-Stokes  Equations  in  Generalized  Curvilinear  Coordinate  System 

Using  a  Cartesian  coordinate  system  to  adequately  describe  flows  with  complex  topology  when  solving  problems  by  finite 
difference  methods  is  difficult  for  two  reasons:  interpolation  procedures  are  needed  to  obtain  boundary  conditions;  there 
are  difficulties  in  describing  the  computational  grid,  since  the  computational  grid  is  not  rectangular.  So,  let  us  move  on  to 
an  arbitrary  curvilinear  coordinate  system.  The  uniform  grid  for  the  difference  approximation  of  the  initial  equations  will 
be  set  in  the  following  coordinate  system: 

T  =  tt  £  =  £{t,x,y,z)t  rj  =  Jj{t,x,y,z )  C  =  £{t,x,y,z) 

The  use  of  generalized  transformation  allows  to  keep  strictly  conservative  form  of  equations.  The  equations  take 

the  form 

AL  +  A*  +  ^  +  ^  =  0 

drj  dfj  drjJ  d(,J 

where 


E  =  ZtU  +  fx(E  -  Ev)  +  (y(F  -  Fv)  +  UG  ~  Gv), 
F  =  r]tU  +  r]x{E  -  Ev)  +  rjy{F  -  Fv)  +  rjz(G  -  Gv), 
G=<{tU  +  CX(F  -  Ev)  +  (y(F  -  Fv)  +  (Z(G  -  Gv). 


The  coefficients  of  the  transformation  matrix  can  be  calculated  by  the  following  formulas: 

r  _  j  ,dy  dz  dz  dy^  ^  _  j.dzdx  dx  dz^  ^  _  j  ,dx  dy  dy  dx^ 

^-^drjdS  drj  dC’ ~  J^dri  dZ  drj  d(J ’  ~  ^drj  drj  d(J’ 

_  ,dz  dy  dy  dz^  _  ,dz  dx  dx  dz.  _  ,dx  dy  dy  dx^ 

ar’  ar’  ~J^  ac 

7  —  jrdz  Qy  dy  dz^  ,  _  .dz  dx  dx  dz^  ,  _  ,dx  dy  dy  dx^ 

Sx  _  ]\j%d^  ~  ^y  -]\)%d^  ~  ĕ^d^J’  z  ~J^dĕ  -  ĕ^ĕnJ,‘ 

where  /  -  is  the  Jacobian  of  the  transformation,  which  is  defined  by  the  formula 

J  d%  dr\  d%  dr\  dr)  dr)  dr\  dr) 


The  use  of  generalized  transformation  allows  to  construct  a  uniform  grid  in  the  form  of  a  unit  cube.  The 
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coefficients  of  the  transformation  matrix  can  be  calculated  either  analytically  (if  the  transformation  is  given  algebraically) 
or  using  difference  formulae  according  to  equations  (if  the  distribution  of  nodes  in  the  physical  domain  of  the  calculation  is 
given).  The  construction  of  a  grid  in  a  physical  domain  adapted  to  real  geometry  is  currently  an  intensively  developed 
problem  of  computational  mathematics.  Suppose  that  the  grid  is  constructed  so  that  near  the  surface  of  the  body,  the 
coordinate  lines  ([xi]  =  const.,  \zeta\  =  const.)  are  oriented  along  the  normal  to  the  surface  [eta\  =  const. 

When  writing  equations,  it  is  assumed  that  the  derivatives  that  are  present  in  the  expressions  for  Ev,  Fv  and  Gv 
are  transformed  according  to  the  rules  of  differentiation  of  complex  functions.  These  components  are  responsible  for  the 
presence  of  viscous  forces.  A  wide  class  is  represented  by  flows  in  which  the  viscosity  forces  significantly  affect  the  flow 
only  in  thin  layers  or  near  the  surface  of  the  streamlined  body.  In  the  case  of  weak  visco-inviscid  interaction  (there  is  no 
interaction  of  the  boundary  or  shear  layer  with  the  incident  jumps  of  compaction  and  separation  of  the  flow).  It  is  possible 
to  solve  the  problem  by  the  classical  approach,  which  consists  of  dividing  the  problem  into  inviscid  and  viscous. 

2.3.  System  of  Euler  Equations 

The  system  of  Euler  equations  in  the  generalized  coordinate  system  is  obtained  from  the  above  system  by  the  limit 
transition.  Viscous  terms  are  assumed  to  be  zero.  The  system  of  Euler  equations  has  the  form 

£.£  +  i.£  +  i.£  +  i.£  =  o 

9t  J  3f ;  drt  J  dt,  ] 

where 

E  =  ftU  +  fxE  +  fyF  +  fzG, 

E  =  r]tU  +  rjxE  +  r]yF  +  r]zG, 

G=<;tU  +  ZXE  +  {yF  +  &G. 

If  we  enter  the  contravariant  components  of  the  velocity  vector 

u  =  C  +  %xu  +  v  A  =  T]t  +  T]XU  +  7] yV  +  T]ZW,  W  =  +  (XU  +  <yV  +  (ZW, 

A  A  A 

the  vectors  E,  F,G  can  be  represented  as 

pu  p v A  pw 

puu-\-  f xp  p  u  v  A  +  TjxP  puw  +  ZXP 

Ĕ  =  [pvu  +  tyP},F  =  {PwA  +  r7yp}?G  =  {pvw  +  (yp}. 
pwu  +  ^yp  pwvA  +  r]yP  pww  +  ^p 

{e  +  p)u  {e  +  p)  v  A  {e  +  p)  w 

2.4.  System  of  Navier-Stokes  Equations  in  the  Thin  Layer  Approximation 

For  the  calculation  of  flows  with  strong  viscous-inviscid  interaction  is  considered  a  single  equation  that  takes  account  of 
both  inviscid  and  viscous  forces.  It  is  assumed  that  the  dissipative  processes  in  the  direction  of  the  normal  to  the  surface  of 
the  body  or  shear  layer  are  much  more  dissipative  processes  in  tangent  directions.  Hence,  in  the  Navier-Stokes  equations, 
when  calculating  dissipative  terms,  it  is  possible  to  neglect  terms  with  derivatives  in  tangent  directions.  As  a  result,  we 
obtain  equations  that  contain  all  the  terms  of  the  Euler  equations,  the  boundary  layer  equations,  and  the  second  derivative 
in  the  normal-to-surface  momentum  equation.  The  resulting  system  of  equations  will  later  be  called  the  Navier-Stokes 
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system  of  equations  in  the  thin  layer  approximation.  This  system  has  the  form 


dr  J  dU  dri  J  dt  J  d7] 

where 


0 

du 

—  +  m2rix 

S  =  i^+m^ 

dw  , 

mi  —  +  m2r]z 
m3 


m!=Vx+Vy+  V  l 


m2 


1  r  du  ,  dv  ,  dw-. 

[Vx^  +  Vy7Z+Vz—l 


r  V  1  dT 

mo  =  mi - —  + 

6  1  Ly-1  Pr  dri 


3  L  dr]  W  d7] 
d  u2+v2+w2 


d7]A 


dr\  d7] 


]  +  m2[r]xu  +  %v  +  rjzw]. 


rj  =  const.  sets  the  boundary  -  the  surface  of  the  body. 

Formally,  the  Navier-Stokes  equations  in  the  "thin"  layer  approximation  are  obtained  by  estimating  the 
derivatives  d/d^  (tangent)  and  d/dlj  (normal)  in  the  viscous  terms  and  excluding  all  derivatives  d/d£ . 


When  solving  problems  of  aerodynamics  in  many  important  cases,  it  is  necessary  to  take  into  account  the  effects  of 
viscosity.  The  separation  of  the  flow  and  the  formation  of  vortex  structures,  the  interaction  of  shock  waves  with  the  boundary 
layer  fundamentally  changes  the  flow  structure,  which,  in  turn,  determines  the  pressure  distribution  and  the  total  aerodynamic 
properties.  In  many  modern  aircrafts,  vortex  structures  are  used  purposefully  to  achieve  certain  aerodynamic  properties  in  a 
wide  range  of  angles  of  attack  and  dellection  of  the  control  of  aerodynamic  surfaces.  For  this  reason,  it  is  often  not  possible  to 
make  a  classical  partition  of  the  problem  into  inviscid  (within  the  Euler  equations)  for  determining  the  pressure  distribution 
and  viscous  (within  the  boundary  layer  equations)  for  estimating  friction.  It  is  necessary  to  use  more  complex  mathematical 
models.  However,  when  solving  three-dimensional  problems  of  extemal  flow  at  high  temperatures,  the  possibilities  of  real 
computer  technology  are  limited  to  the  condensation  of  nodes  only  in  one  coordinate  direction.  The  greatest  gradients  due  to 
viscosity  are  realized  in  a  thin  viscous  layer  near  the  surface  of  the  body.  Accordingly,  the  thickening  of  the  nodes  is  carried 
out  in  the  direction  of  the  normal  to  the  surface  of  the  body,  which  makes  it  possible  to  calculate  the  terms  of  the  boundary 
layer  equations.  For  this  reason,  from  the  point  of  view  of  real  computational  possibilities  in  calculating  the  dissipative  term  of 
the  Navier-Stokes  equations,  it  seems  reasonable  to  limit  the  calculation  of  the  terms  included  in  the  boundary  layer 
equations,  namely  only  the  second  derivatives  of  the  normal  to  the  surface  of  the  body. 

This  model  is  called  the  Navier-Stokes  equations  in  the  thin  layer  approximation.  This  model  should  be 
distinguished  from  the  thin  shock  layer  model,  which  is  used  to  solve  hypersonic  aerodynamics  problems  near  blunted 
bodies. 


The  system  of  Navier-Stokes  equations  in  the  thin  layer  approximation  contains  all  the  terms  of  the  Euler  and 
boundary  layer  equations,  which  allows  to  take  into  account  both  the  inviscid  and  the  main  viscous  effects,  and  as  a  result, 
the  visco-inviscid  interaction  within  the  framework  of  unified  equations.  The  use  of  this  model  is  mainly  due  to  the 
resources  of  existing  computer  calculation. 

2.5.  Turbulence  Models 

The  calculation  of  turbulent  flows  is  carried  out  in  the  framework  of  averaged  Navier-Stokes  equations  with  different 
closure  models:  algebraic  model  from  (Baldwin  and  Lomax,  1978),  model  from  (Tulapurkara,  1997;  Belov  and  Isaev, 
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2001),  model  of  time-averaged  kinetic  energy  of  turbulent  motion  and  turbulent  yiscosity,  LARGE-scale  Les  yortices, 
direct  numerical  simulation  (DNS). 

2.6.  Numerical  Method 

We  brietly  describe  an  explicit  second-order  difference  scheme  of  predictor-corrector  approximation  (see  MacCormack, 
1969).  We  introduce  a  grid  with  step  sizes  h\,  h2  and  /i3?  such  as  that  of  =  mhh  r\k  -  kh2  and  Q  =  lh2.  Integration  in  time  to 
achieye  a  stationary  state  is  carried  out  on  the  basis  of  a  two-step  scheme,  which  is  formally  represented  as 


u!l!,,,  =  u"  -  f  (e:  -  e;_,  )-  f  (f;  -  f;_,  )-  f  (g;  -  g;_,  ) 


u;:;,=||u"+u(,|-|-(E;),-E«)--^(F«,-F,(i|)-^(G!;,-G!,>)j 


Here,  the  superscript  (1)  denotes  that  the  function  values  are  calculated  for  intermediate  solution  values  U(1).  In 
addition  to  the  written  difference  relations,  the  next  final  step  is  the  “smoothing”  of  the  solution,  often  necessary  in  areas 
with  shock  waves. 

In  fact,  convective  terms  are  represented  by  one-sided  differences  with  a  change  of  direction,  and  dissipative 
terms  are  represented  by  symmetric  differences. 

When  using  explicit  difference  schemes,  the  time  step  size  is  limited  by  the  stability  condition  of  the  calculation; 
we  use  the  method  of  local  linearization. 

In  the  case  of  compressible  gas  flows,  the  Euler  equation  system  is  written  as 

au  Aau  Dau  _au  n 

- +  A —  +  B —  +  C —  =  0 

dT  J  dg  dTJ  dC 


Eigenvalues  of  a  matrix  A,  B  and  C: 


°lX3=B  =  TlxU+7lyV  +  T]zW 


(71C2,3  =C  =  CxU  +  CyV  +  CzW 


rji 

H 

<5=c±  U(£+cj+£). 


The  restriction  on  the  time  step  is  represented  in  the  form 
At  <  Ku — At  <  Ku — -p-g-,  At  <  Ku — 

max|crA|  max|<T°|  max|cr^| 


where  Ku  is  the  Courant  number. 


www.tivrc.ors 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


642 


Vassiliy  Kryuchkov,  Vitaly  Savinkin,  Andrey  Solodovnik  &  Lyudmila  Kryuchkova 


Viscous  dissipative  effects  have  a  great  inAuence  on  the  flow  near  the  surface  of  the  body.  In  addition,  actually 
used  computational  grids  due  to  limited  computational  resources  do  not  provide  the  possibility  of  sufficient  thickening  of 
nodes  at  a  distance  from  the  body.  For  these  reasons,  the  dissipative  term  can  be  taken  into  account,  and  therefore  the 

calculation  can  be  carried  out  only  in  the  region,  where  0  <  7j  <  7j  and  7j  <  1  (rj  =  1  -  the  outer  boundary),  and  Tj 
determined  either  experimentally,  or  a  condition  must  be  set  for  its  determination.  In  the  region  7j  <7j  <  1  the  dissipative 
term  is  small  and  can  be  put  equal  to  zero. 


The  introduction  of  the  dissipative  term  changes  the  stability  condition.  The  condition  of  the  restriction  on  the 
step,  taking  into  account  the  dissipative  term  has  the  form: 


At  <  Ku 


Ar) 


|<TB|+2|crv|/ar7 


,  meav  =  (ril  +  Vy  +  Vz) 


7M 

pRe  Pr’ 


Note  that  accounting  for  viscosity  in  the  complete  or  simplified  Navier-Stokes  equations  requires  a  strong 
thickening  of  nodes  in  viscous  thin  layers,  which  leads  to  a  significant  reduction  in  the  integration  step  over  time,  but  this 
limitation  is  mainly  due  to  the  rate  of  propagation  of  perturbations,  not  viscosity.  The  type  of  the  step  restriction  taking  into 
account  the  dissipative  term  reveals  that  with  a  fixed  grid,  an  increase  in  Re  leads  to  the  possibility  of  increasing  the 
integration  step.  If  the  Re  decreases  (i.e.,  the  viscosity  of  the  flow  increases),  the  possible  integration  step  decreases.  In 
reality,  the  grid  step  is  related  to  the  Reynolds  number.  The  greater  the  Re,  the  greater  must  be  the  thickening  of  the  nodes 
to  describe  thin  viscous  layers,  and  the  reduction  of  Re  allows  for  the  use  of  grids  with  less  thickening.  In  fact,  the  design 
grid  should  be  built  with  Re  in  mind. 


To  perform  the  calculation,  the  Courant  number  is  usually  assumed  to  be  Ku  from  0.5  to  0.9.  Due  to  the  need  to 
construct  a  substantially  uneven  grid,  the  integration  step  in  time  is  determined  in  the  region  of  the  greatest  condensation  of 
nodes.  Establishing  a  single  integration  step  is  extremely  slow.  A  local  integration  step  is  used  to  speed  up  the  installation: 
the  integration  step  at  each  node  is  selected  locally  based  on  local  conditions.  This  ensures  that  the  distribution  of  the 
disturbance  with  a  speed  of  one  knot  for  one  integration  step  with  Ku  =  1.0.  Using  a  local  integration  step  leads  to  the  need 
to  reduce  the  number  of  chimes  from  0.1  to  0.25.  The  use  of  the  local  integration  step  does  not  allow  to  consider  non- 
stationary  problems,  but  if  only  the  steady-state  flow  is  of  interest,  then  the  local  step  allows  to  significantly  reduce  the 
total  number  of  integration  steps  before  the  solution  is  established. 


The  explicit  scheme  described  has  second-order  accuracy  in  time  and  spatial  variables  on  smooth  solutions. 
Changing  the  pattem  or  increasing  the  number  of  iterative  steps  in  time  can  increase  the  order  of  approximation  in  space  or 
time. 


The  scheme  of  the  second  order  in  accuracy  in  time  and  the  fourth  in  space  for  a  one-dimensional  unsteady 
equation  has  the  form 

Um  =  Un  -  iU7(££+l  -  Em)  -  (Em+ 2  ~  Em+l)l 

Um  =  \{Un  +  t/(1)  -  [7 (41}  -  Em-l)  ~  (Em-1  ~  Em-2)B- 

Increasing  the  order  of  approximation  in  time  is  achieved  by  introducing  additional  time  steps  (Runge-Kutta  type 
schemes): 
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jj^  —  jjn  _  i_L  (J7n  _  p n  \ 
um  ~u  2  h± 


643 


U^  =  un---(E (1) 

um  u  9  ^  Vnm+ 1  nm  J 


2  ht 


rr(^)  _  jjn  _  ( /7(2)  _  /7(2)  ^ 

Um  ~u  2hx^m  ^m-lJ’ 


Um1  =  Un-  ~{(En  -  En_J  +  2 (£■«,  -  E%>)  +  2 (42)  -  42^)  +  (43^  -  4°)}. 


The  use  of  explicit  schemes  allows  a  simple  implementation  of  a  numerical  algorithm  for  solving  differential 
equations  on  multiprocessor  computer  technology. 


2.7.  The  Process  of  Designing  a  Model  and  Creating  a  Computational  Grid 

To  solve  the  flow  problem,  it  is  necessary  to  describe  the  geometry  of  the  object  of  study.  The  construction  of  the  geometry 
depends  on  the  initial  information  and  the  objectives  of  the  study.  At  all  stages  of  design,  the  choice  of  a  real  form  is  one  of 
the  central  tasks.  There  are  a  number  of  computer  applications  that  allow  you  to  automate  geometric  design.  An  important 
feature  of  these  systems  is  the  ability  to  directly  observe  the  results  of  the  work  on  the  computer  screen.  In  the  construction 
of  complex  geometry,  various  methods  are  used:  "solid-state"  method,  analytical  description  and  spline  theory.  Since  each 
system  has  its  own  goals  and  uses  original  methods  of  describing  geometry,  the  transfer  of  information  from  one  system  to 
another  is  difficult.  For  data  exchange  between  automation  systems,  special  data  formats  have  been  developed,  focusing  on 
which  and  having  the  means  of  their  analysis,  it  is  possible  to  obtain  the  necessary  geometry  data  (for  example,  data  on  the 
surface  of  a  streamlined  body,  perpendicular  to  its  surface). 

Solid-state  geometric  modeling  system  consists  of  a  set  of  tools  that  allow  you  to  interactively  create  from  simple 
geometric  shapes  (plane,  cone,  cylinder,  ball,  ellipsoid,  etc.).  The  geometry  of  the  real  form  using  the  operations  of  union, 
intersection  and  subtraction.  The  division  of  geometry  into  bodies  of  simple  form  and  the  subsequent  construction  of 
geometry  consists  of  determining  the  geometric  parameters  of  simple  bodies,  determining  their  position  in  space  and  then 
combining.  Such  applications  of  programs  for  work  with  geometrical  models  allow  to  design  and  create  objects 
independently,  using  the  saved-up  experience  from  a  database,  to  make  geometry  of  aerodynamic  profiles,  using  sets  of 
various  data  (for  example,  radius  of  curves,  angles  of  inclination  of  various  axes,  radii  of  a  fuselage,  angle  of  attack  of  a 
wing,  degree  of  placement  of  a  tail  part,  radii  of  sides,  etc.)  from  volume  primitives.  Computer-aided  design  technology 
based  on  the  features  of  models  (figure  1),  makes  it  possible  to  increase  the  speed  of  design,  improve  the  accuracy  of  the 
values  of  the  studied  processes  as  well  as  improve  the  aerodynamic  characteristics  of  the  designed  airframe,  which  will 
largely  reduce  the  complexity,  material  consumption,  error  in  the  design  and  facilitate  the  creation  of  a  new  aircraft  with  an 
optimized  shape. 

Also,  in  the  development  of  design,  offices  and  research  institutes  used  various  methods  of  creating  models  based 
on  three-dimensional  scanning  of  real  samples  (figure  2).  However,  such  models  have  a  large  number  of  known 
disadvantages  -  time-consuming  surface  preparation,  the  presence  of  insensitive  areas  for  the  scanner,  the  complexity  of 
processing  a  large  number  of  roughness  on  the  surface,  etc. 
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Figure  1:  UAY  Model  Built  using  Solid-State  Design. 


Figure  2:  UAV  Model,  Built  using  Three-Dimensional  Scanning 
(Reyealed  Only  the  Fuselage). 


As  can  be  seen  from  the  figures,  the  most  acceptable  model  for  calculations  is  precisely  built  because  of  its 
symmetry  and  better  streamlining. 


When  solving  flow  problems  numerically,  it  is  important  to  construct  not  only  the  geometry  of  the  body,  but  also 
calculation  of  grid  near  it.  In  order  to  include  grid  construction  in  the  design  system,  you  need  to  “correctly”  construct  the 
geometry  using  the  computer-aided  design  system.  Building  a  grid  around  arbitrary  bodies  is  a  time-consuming  process 
that  depends  on  many  factors  arising  in  the  course  of  the  project.  Quite  often  such  surfaces  are  given  from  the  point  of  view 
of  parameterization  unsuccessfully.  Therefore,  not  only  must  the  construction  of  the  grid  be  a  sustainable  process,  but  the 
construction  method  should  not  depend  on  the  method  of  parameterization.  For  optimization  problems,  design  parameters 
must  be  selected,  so  that  they  are  sensitive  to  grid  changes. 


When  solving  problems  of  external  flow  by  methods  of  computational  aerodynamics,  the  calculated  grid  must  be 
connected  to  the  surface  of  the  streamlined  body,  which  is  due  to  the  nature  of  the  flow  parallel  to  the  surface  of  the  body. 
In  addition,  when  taking  into  account  the  viscosity  to  describe  the  solution,  it  is  necessary  to  significantly  thicken  the  nodes 
to  the  surface  of  the  body  to  describe  the  profile  of  gas-dynamic  functions  in  the  thin  layer  near  the  body.  If  these  features 
are  not  taken  into  account,  it  seems  impossible  to  build  a  solution  that  is  adequate  to  reality.  However,  in  the  external  flow 
region,  it  is  advisable  to  use  a  more  regular  computational  grid  that  is  not  connected  to  the  surface  of  the  streamlined  body. 
We  develop  methods  of  constructing  meshes  to  allow  to  some  extent  to  smooth  out  the  effects  of  fractures  of  the  body 
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surface  to  the  mesh  when  removed  from  it,  but  the  problems  associated  primarily  with  the  complexity  of  the  topology  of 
the  considered  task,  leads  to  the  need  to  build  a  multiblock  mesh,  sew  on  some  borders.  This  complicates  the  algorithms  for 
solying  problems.  It  is  proposed  to  build  computational  grids  near  a  streamlined  body  (or  several  bodies)  and  in  an  external 
flow  independently.  In  addition,  the  grid  near  the  streamlined  body  is  built  only  in  some  neighborhoods  of  the  body,  which 
greatly  simplifies  its  construction. 

The  calculated  main  grid  (figure  3)  is  a  grid  with  rectangular  cells  and  is  designed  to  calculate  the  external 
inviscid  flow. 


Figure  3:  Construction  of  the  Calculated  Grid  around  the  Body  Under  Study. 


Thus,  during  the  simulation,  we  choose  a  high-speed  coordinate  system,  which  is  tied  to  the  airspeed  of  the 
aircraft.  In  this  case,  the  speed  of  the  vessel  is  defined  as  a  vector  consisting  of  components  Vx,  Vy,  Vz  : 

\  =  Vxi+Vyj+Vzk  ,  and  it  gives  us  direct  information  about  the  strength  of  the  frontal  resistance  and  lifting  force. 

For  the  calculation,  we  will  take  the  following  parameters  of  the  external  environment  according  to  the  model  of 
the  standard  atmosphere  (Altitude  4000  m)  (Anderson,  2016): 

Medium  type:  Air 

Flight  height:  4000  m 

Atmospheric  pressure  at  this  altitude:  61300  Pa 
Temperature:  -100°C 
Center  of  mass  speed:  20  m/s 
Angle  of  attack:  [alpha]  =  6  degrees 

The  following  projects  were  selected  as  the  models  under  consideration  (figures  4-9): 

•  The  basic  model  without  changes. 

•  A  basic  airframe  with  a  single  spherical  dome  fairing  mounted  at  the  front  of  the  fuselage. 

•  Basic  airframe  with  a  single  spherical  dome  fairing  mounted  in  the  middle  of  the  fuselage. 

•  Basic  airframe  with  two  spherical  dome  fairings  of  different  diameters,  mounted  on  the  front  and  middle  of  the 
fuselage. 
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•  Base  glider  with  combined  domes  of  different  diameters. 

•  The  basic  airframe  with  an  elongated  fairing  of  elongated  shapes. 

As  output  parameters,  we  will  compare  changes  in  resistance  forces  from  air  flows,  lift,  changes  in  flow 
temperature  (minimum,  average  and  maximum)  as  well  as  moments  of  forces  acting  on  the  airframe. 


Figure  4:  UAY  Profile  of  Project  1. 


Figure  5:  UAY  Profile  of  Project  2. 


Figure  6:  UAY  Profile  of  Project  3. 


Figure  7:  UAY  Profile  of  Project  4. 


Figure  8:  UAV  Profile  of  Project  5.  Figure  9:  UAV  Profile  of  Project  6. 

After  that,  a  mathematical  calculation  was  carried  out  according  to  the  mathematical  model  described  above. 
Based  on  the  results  of  calculations  to  visualize  changes  in  air  flows  and  pressure  redistribution,  pressure  distribution 
patterns  were  constructed  in  the  plane  passing  through  the  central  axis  of  symmetry  of  the  airframe  (figures  10-15),  which 
allow  to  see  the  result  of  changes  in  aerodynamic  shape. 
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Figure  10:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  1. 


Figure  11:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  2. 


Figure  12:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  3. 
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Figure  13:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  4. 


Figure  14:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  5 


Figure  15:  Pressure  Isobars  in  the  Symmetry  Plane  of  Project  6. 


3.  RESEARCH  RESULT 

A  detailed  analysis  of  the  pressure  changes  revealed  that,  indeed,  changes  in  the  upper  fuselage  did  not  have  a  significant 
change  in  the  pressure  distribution  in  the  lower  airframe.  However,  there  are  significant  changes  in  the  upper  hemisphere 
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of  the  surrounding  space.  The  reason  for  this  change  lies  in  the  change  in  the  current  lines  Aowing  around  the  airframe.  To 
do  this,  we  also  constructed  the  trajectories  of  the  current  lines  of  the  air  particles  depending  on  their  speed  and  location. 


Figure  16:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  1. 


Figure  17:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  2. 


Figure  18:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  3. 


Figure  19:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  4. 
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Figure  20:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  5 


Figure  21:  Yelocity  Distribution  of  Air  Current  Lines  of  Project  6 


For  global  goals,  data  tables  were  created  for  the  values  of  interest  for  each  of  the  projects  and  as  a  result,  a 
common  comparatiye  table  1  was  formed,  which  is  given  below. 


Table  1:  Comparatiye  Characteristics  of  Parameter  Changes 


Project  1 

Project  2 

Project  3 

Project  4 

Project  5 

Project  6 

The  Lifting  Force  Fl,  N 

83,37 

100,22 

94,16 

89,28 

94,23 

85,97 

Resistance  force  Fr,  N 

6,26 

0,17 

6,38 

7,25 

4,29 

6,25 

Pitch  Moment  Mx,  N*m 

10,25 

11,01 

11,02 

10,14 

11,10 

10,29 

The  air  temperature,  K 

263,06 

263,09 

263,07 

263,08 

263,08 

263,09 

From  the  data  of  the  table,  it  can  be  concluded  that  the  parameters  of  air  temperature  change  when  changing  the 
aerodynamic  shape  cannot  be  taken  into  account,  because  the  change  occurs  only  by  hundredths  of  a  degree. 

Analyzing  the  table  dependence,  it  can  be  concluded  that  the  latest  modification,  presented  as  project  6, 
practically  does  not  make  changes  in  the  moments  of  pitch  forces  and  drag  forces  Fr  and  lifting  F1  forces.  In  some  cases, 
the  drag  force,  as  expected,  increased,  which  can  be  observed  in  the  visualization  of  air  flows  around  the  airframe. 
Interestingly,  behaves  lift  Fl,  which  slightly  increased  by  20%  after  the  first  modification,  but  then  gradually  decreased  to 
86  N.  But  the  most  significant  results  of  the  study  are  changes  in  the  moments  of  forces  -  Mx,  My,  Mz  -  acting  on  the 
projected  airframe.  The  pitch  moment  Mx  about  a  horizontal  axis  perpendicular  to  the  direction  of  movement  of  the 
airframe  was  changed  (and  by  how  much?)  by  increasing  the  drag  Fr,  but  it  can  be  compensated  by  dellecting  the  Aaperons 
within  up  to  5  degrees  before  their  significant  contribution  to  the  drag  force  occurs,  which  will  entail  some  change  in 
aerodynamic  shape.  Changing  the  design  geometry  and  shape  will  ultimately  affect  the  duration  of  the  Hight.  These 
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Thus,  we  have  proven  that  at  an  angle  [alpha]  =  6  degrees,  the  characteristic  of  the  horizontal  flight  of  most 
aircrafts,  the  latest  modification  with  a  single-lamp  aerodynamic  profile  is  the  most  optimal  for  the  implementation  of  our 
project. 

Let  us  explore  the  latest  modification  and  see  how  the  lift  F1  and  drag  force  Fr  behave  depending  on  the  angle  of 
attack.  We  will  vary  it  from  0  to  10  degrees. 

The  constructed  graphs  of  dependencies  are  presented  in  figures  22-24. 


Dependence  of  the  Lift  Force  on  the  Angle  of  Attack 
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Figure  22:  Graph  of  the  Dependence  of  the  Lift  Force  (in  N)  on  the  Angle  of  Attack  (in  u). 
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Figure  23:  Graph  of  the  Dependence  of  the  Drag  Force  (in  N)  on  the  Angle  of  Attack  (in  °). 


Figure  24:  Graph  of  the  Dependence  of  the  Pitch  Moment  (in  N*m)  on  the  Angle  of  Attack  (in  °) 
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4.  CONCLUSIONS 

As  can  be  seen  from  the  graphs,  with  an  increase  in  the  angle  of  attack  [alpha]  from  0  to  10  degrees,  the  lift  FL  and  the 
moment  of  pitch  Mx  increase  practically,  as  predicted  in  the  work  (Fundamentals  of  aerodynamics  and  Aight  dynamics, 
2010).  It  is  also  concluded  that  the  aerodynamic  profile  of  the  proposed  model  6  (figure  9)  is  not  symmetrical,  because  the 
graph  of  the  lift  dependence  (figure  22)  does  not  start  from  the  origin  (not  from  the  point  (0,0)).  The  pitch  moment  (figure 
24)  also  grows  from  3  to  16  N*m  with  increasing  angle  of  attack  from  1  to  10  degrees,  almost  linearly.  This  process  is 
explained  by  the  fact  that  the  difference  in  the  total  forces  applied  to  the  nose  and  tail  parts  of  the  airframe  increases  due  to 
the  different  surface  area.  Special  attention  should  be  paid  to  the  result  that  we  see  in  figure  23,  where  there  is  a  clear 
minimum  resistance  force,  Fr  =  2.14  N,  at  the  angle  of  attack  [alpha]  =  5  degrees,  in  contrast  to  the  standard  values  of  the 
angle  of  60  degrees.  In  this  regard,  it  can  be  concluded  that  the  addition  of  the  upper  fairing  led  to  a  change  in  the  ahTrame 
fairing,  which  must  be  taken  into  account  when  changing  the  autopilot  programs. 
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